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METHODS OF OPTIMISZING STOCHASTIC PROCESSING 
PARAMETERS IN CROP HARVESTING MACHINES 

CROSS-REFERENCE TO RELATED PATENT APPLICATIONS 

fOOOn This Patent Application claims priority under 35 U.S.C. 371 to PCT 
International Application Number PCT/EP03/08272. which was filed on July 25, 
2003. and which has Convention Priority based on Great Britain Application GB 
02 1 7297. 1 . filed July 26. 2002> 

FIELD OF THE INVENTION 

100021 T his invention relates to methods for establishing parameters for 
modelling the behaviour of crop harvesting machines such as but not limited to 
combine harvesters. In particular it relates to a method for optimizsing stochastic 
model parameters in models for a crop processing operation in such harvesting 
machines. 

The method is particularly useful in the context of "precision farming" techniques, 
in which combine harvesters and other harvesting machines have played a 
significant role, as well for the development as for the application of these 
techniques. 

BACKGROUND OF THE INVENTION 
[00031 " Precision farming" is a broad term that embraces practices such as: 



30 - field yield mapping aimed at accurately mapping the crop-producing 
productivity of a field so that: 

• seeds and agro-chemicals can be economically and correctly applied 
without over- and under-dosing the plants; 

1 



20 



25 



• harvesting and other agricultural machinery can be adjusted to take 
account of varying crop conditions from place to place in a field; 

- modelling, of process conditions, within agricultural machines, and developing 
5 control philosophies aimed at: 

• improving machine efficiency or workrate; 

• reducing fuel consumption; and 

I • improving the quality of work carried out by agricultural machines; and 

10 I - providing alerts and reports of abnormal conditions in crops-. 

I 100041 Sometimes the overall improvement in farming efficiency resulting 
from an individual precision farming practice might be only a few percent; but 
cumulatively such efforts have dramatically improved productivity in mechanised 
15 farming over recent years. 

I rOOOSI T he significance of harvesters, such as combine harvesters, to precision 
farming derives principally from the following factors: 

20 (i) A harvester is able to assess the output of a farming operation for example in 
terms of crop yield. A field map of such data is invaluable in improving 
I farming efficiency in fixture crop growing and harvesting seasons: and 

(ii) Harvesters are complicated machines that provide numerous sites for the 
25 location of transducers whose fiinction is to gather data on the harvesting 

processes and the crops passing through the machine. 

I 10006] A s noted, the combine harvester has been particularly usefiil in 
producing maps indicating the expected crop yield at different locations in a field. 
30 Farmers can use such maps (that are readily stored in digital form in a computer 
memory) to control in an accurate way the processes forming part of a crop 
growing season, so that the yield of the field is maximizsed. 
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[00071 Prior art techniques for yield mapping, however, are limited primarily 
because they concentrate on the quantity of the useable part of the crop that is 
conveyed to e.g. the clean grain tank in a combine harvester. 



fOOOS I A lthough measurements of e.g. the mass flowrate of clean grain to the 
clean grain tank may readily be compensated for some variables such as grain 
moisture content and grain type, some difficulties remain. 



j [00091 Significant among these is the fact that mass flowrate measurements of 
10 crop yield generally take no account of crop losses arising from e.g. incomplete or 
faulty threshing of ears in the threshing drum or another part of the harvester 
where grain separation occurs. Where the machine load exceeds the threshing, 
separation and/or cleaning capacity of the machine a portion of the harvested grain 
will not be separated from the straw and chaff and be deposited therewith on the 
15 field behind the machine. 



I rOOlOl I t is known, for example, that the extent to which (or the ease with 
which) ears are threshed in the threshing drum of a combine harvester is strongly 
dependent on the feedrate of crop into the harvesting machine, when such factors 
20 as grain and straw moisture, crop variety and straw length are kept as constants. 
Thus for higher feedrates relatively less grain is separated in the first concave 
threshing drum section than in the case of low feedrates. More grain has to be 
separated in the further stages of the harvester. 

25 [00111 Consequently for such high feedrates a lesser proportion of the crop is 
therefore likely to reach the clean grain tank, with the result that a yield 
measurement taken at such a location may be inaccurate. Furthermore such a 
measurement takes no account of the extent to which grains become damaged or 
lost within the combine harvester. 
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[00121 I n reality a great number of variables influences the extent to which the 
threshing and separating sections are able to separate grain from other plant matter 
such as chaff and straw. Such variables include, but are not limited to: 
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- the nature of the soil in which the crop grows; 

- settings of various adjustable components of the harvesting machine, e.g. the 
height of the header bar in a combine harvester, which directly influences the 
straw to grain ratio; 

5 - the slope of the field in which the machine operates; 

- the moisture content of the crop; 

- the crop type; 

- the forward speed of the harvesting machine; 

- the presence of weed patches; 

10 - the state (wear) of the machine elements; 

- the type of installed machine elements, e.g. the type of rasp bars; 
and so on. 

room A ccording to a first aspect of the invention there is provided a method 
of substantially continuously optimizsing a stochastic parameter 9 that 
characterizses the instantaneously prevailing readiness with which crop is 
processed in a harvesting machine, including the step of recursively calculating 
the optimizsed parameter value in accordance with the following algorithm: 

m-f{kt-\)MtM'-m - (a) 

wherein: 

&(t) is the optimizsed stochastic parameter value at time t; and 
£{t, i9(/)) is an error prediction function. 

[00141 Such a method is highly suited to the continuous optimizsation of the 
25 highly stochastic parameter i9 that, when applied to the threshing and separation 
process a combine harvester, may fairly be termed a "threshability" parameter, i.e. 
an indication of the extent to which the harvesting machine is capable of threshing 
the crop at time t, 

30 I [00151 Such a parameter is useable in various ways, as discussed hereinbelow. 

I f00161 T he algorithm generally may take the form of: 
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j r00171 T he method of the broad aspect of the invention can readily be carried 
out using a suitably programmed computer carried by or forming part of the 
5 harvesting machine. 

I f 00 1 81 Preferably the algorithm (A) has the form-: 

s{t) = s{t - 1) + r{t)r-' (/ Vt. - i)M^ - 0) 

wherein 

10 y{t) is a gain term; 

r{t) is a scalar approximation of a Hessian K"(i9) in which Fis indicative 
I which K is a quadratic error criterion; m4 

y/if^S) - ^y^^^^ J in which y{ty&) is an estimation of a value indicative of 
(19 

the effectiveness of said crop processing in said harvesting machine said 
15 estimation being based on stochastic parameter 9 ; and 

e{t,9{t'-\)) is the difference between the actual effectiveness vdAxxQ y(t) 

and the estimated value y{f,9) based on the previously optimizsed parameter 

^(/-l). 
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[00191 P referably the algorithm (A) includes an estimation of r{t) that is 
weighted to reduce the influence, on the optimizsed parameter values 9 , of past 



measurements. 

I [00201 T his aspect of the method renders the parameter optimizsation more 
25 realistic and robust for a wide range of working conditions. 

I [002 1 1 The parameter 9 may be usable in a model for the relation between a 

value u(t) indicative of the feedrate of crop into the harvesting machine and a 

value y(t) indicative of the effectiveness of an operation processing said crop in 

30 said harvesting machine. The estimated value y{ty9) is then an estimation of the 

effectiveness obtained by the application of said model to the feedrate values u(t). 
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r00221 I n this manner the model can be updated continuously in order to meet 
any changes to the process caused by a wide range of changing conditions, e.g. 
varying crop properties such as ripeness or moisture or changes in inclination of 
the machine. 

[00231 A dvantageously, the model may comprise an exponential function. 
Such form provides some computational advantages for the optimizsation of 3 . 

10 [ 00241 T he effectiveness may take the form of a value indicative of crop flow, 
e.g. crop losses at the end of the separation or the cleaning section. It may also 
comprise the crop flow in a return system. 

I r00251 A ccording to a second aspect of the invention there is provided a 
15 method of operating a harvesting machine comprising the steps of: 

(i) ^substantially continuously optimizing a stochastic parameter that 

characterizses the instantaneously prevailing readiness with which the harvesting 
machine processes crop: and 

(ii) substantially continuously adjusting a performance variable of the 

20 harvesting machine in dependence on the instantaneous, optimizsed value of the 

said parameter in order to optimizse the load of the harvesting machine so as to 
keep a value indicative of the effectiveness of said harvesting machine below a 
predetermined value. 



25 I 100261 Such effectiveness value may comprise the losses of useable crop parts 
such as separation or cleaning sieve losses, or a proportion of damaged useable 
crop parts, e.g. broken grain kernels, or a proportion of unwanted material in the 
useable crop parts, e.g. chaff and straw and particles in the clean grain. 



30 r00271 O ptimizsing the machine load may comprise optimizsing the feedrate 
of crop into the harvesting machine, e.g. by adapting the travel speed of the 
harvesting machine. 
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[00281 C onveniently the step («)-of adjusting a performance variable of the 
harvesting machine occurs in dependence on the output of an inverted form of an 
effectiveness estimation function: 

j)(/,5) = exp(5M(/))-L - (B) 

Herein u(t) may be the measured feedrate and y{t,3) the grain losses. 

[00291 A ccording to a third aspect of the invention there is provided a method 
of mapping one or more field lots for variations in a stochastic parameter that 
characterizses the instantaneously prevailing readiness with which a harvesting 
machine processes crop, the method comprising the steps of: 

^ — operating a harvesting machine to harvest crop in a said field lot; 

(H)-simuItaneously measuring the machine load and the machine 
effectiveness and determining the position of the machine in the field lot; 

(Ui^storing data indicative of the position of the harvesting machine at 

time t; 

(iv) -using the measured machine load and machine effectiveness data in an 
optimizsation of tbe-said parameter; and 

(v) — mapping the optimizsed parameter values obtained from the using the 
measured machine load and machine effectiveness data s tep (f¥)-so as to produce 
a parameter map of the field lot. 

[00301 A ccording to a fourth aspect of the invention there is provided a 
method of operating a harvesting machine comprising the steps of: 

(i) substantially continuously optimising a stochastic parameter that 

characterizses the instantaneously prevailing readiness with which the harvesting 
machine separates useable crop parts from other plant matter: and 

sending a display signal, that is indicative of the instantaneous 

parameter value, to a display device. 



[00311 P referably the display signal indicates an abnormal parameter value. 



[00321 P referably in each of the second, third and fourth aspects of the 
invention the optimizeation step is in accordance with the first aspect of the 
invention. Thus the method of the first aspect of the invention is highly versatile 
5 in its application. 

[003 3 1 Conveniently in each of the 2"** to 4^*^ aspects of the invention, when 
the parameter optimizsation is according to the first aspect of the invention, the 
said selected part of the harvesting machine is selected from: 
10 the separation section, e.g. the straw walkers or a rotary separator; 

the sieve; 

the return flow system; of 
the cleaning section: or 
the grain elevator; 
15 of a combine harvester. 



BRIEF DESCRIPTION OF THE DRAWINGS 



j [00341 T he invention will now be described in more detail, by way of non- 
20 limiting example. In the following description reference is made to the 
accompanying drawing figures in which: 

[00351 Figure 1 is a schematic representation of the threshing, cleaning and 
straw walker sections of a combine harvester showing the possible locations for 
25 grain loss transducing devices; 



[00361 F igure 2 is a graphical representation of the highly stochastic nature of 
a processability, in particular threshability parameter i9, by reference to field 
slope (gradient) values; 

[00371 F igure 3 shows a graphical form the variability of the optimizsed 
parameter value from place to place in a field undergoing harvesting. The Figure 
3 graph also includes, for comparison purposes, plots of parameter values 
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calculated offline using a Least Squares technique also used to generate the plots 
of Figure 2; and 

[00381 F igure 4 is a comparison between a parameter value map and the slope 
5 of the field to which it relates. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

I [00391 T he terms "grain", "straw", and "tailings" are used principally 
10 throughout this specification for convenience and it should be understood that 
I these terms are not intended to be limiting. Thus- "grain" refers to that part of the 
crop which is threshed and separated from the discardable part of the crop material 
which is referred to as *'straw"._ Incompletely threshed ears are referred to as 
"tailings". _Also, the terms "forward", "rearward", "upper", "lower", "left", 
15 "right", etc. when used in connection with the combine harvester and/or 
components thereof, are determined with reference to the combine harvester in its 
normal operational condition and may refer either to the direction of forward 
operative travel of the combine harvester or to the direction of normal material 
flow through components thereof. These terms should not be construed as 
20 limiting. 

[00401 F igure 1 shows the position of respective crop flow or loss sensors 10, 
11, 12 in the separation (13), cleaning (14) and return flow (15) sections of a 
typical combine harvester. The illustrated separation section comprises a set of 
25 straw walkers 13, but the invention may also be used with harvesters comprising a 
rotary separator. 

I [00411 A nv of the positions 10, 1 1 , 12 may readily be used within the scope of 
the invention as locations for grain loss sensors. The following notes apply to the 
30 respective sensors: 

I [00421 Sensor 10: 
Separation losses are generally measured by means of impact detectors. Impacts 
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of kernels are separated from impacts of straw ends and counted each second. By 
this, a measure is found for the amount of separation losses. Since the distribution 
of free grain kernels along the depth of the straw layer varies in time, only a rough 
indication is normally obtained. 

5 

[00431 Sensor 11: 

One similar impact sensor may be installed at the end of the cleaning section. 
Again only a rough measure is normally found for the sieve losses. 

r00441 Sensor 12: 

Tailings fall behind the lower sieve into the cross auger of the return flow system. 
An impeller is installed at the end of the fast rotating auger to spread the return 
flow. The sensor preferably operates measuring the phase shift between a fixed 
impeller and a flexible one as described in EP-A-0 463 240. Higher material 
flows imply larger phase shifts. 

[00451 I n addition, the harvester preferably is equipped with precision farming 
tools to measure the actual cutting width, ground speed, mass flow, feedrate and 
position within the field. 

20 

I [00461 The outputs of such sensors may be employed in a technique, involving 
the use of an offline model for the crop processing operation. Establishing such 
model involves the determination of one or more parameters, e.g. a threshability 
parameter 3 , 

25 

[00471 T he exponential function offers good opportunities to make a 
prediction y{t,S) of actual separation losses y(t) (T/h) based on an online 
measurement of feedrate u(t) (T/h). The feedrate may be measured at the inlet of 
the harvester, e.g. by monitoring the volume of crop material entering the straw 
30 . elevator, or by measuring the force or torque needed to convey the crop mass into 
the machine. 
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[00481 A dvantageously, the exponential function may take the form of: 

10 



j)0,^») = exp(5«(0)-l. 



Equation (1) illustrates that no separation losses occurs for a zero feedrate and the 
5 feedrate-loss relation is fully determined by parameter The time delay 
between the feedrate signal u(t) and y(t) is assumed to be fixed and an optimal 
fixed time shift (typically lis) is installed as a compensation. For online 
measurements, separation losses and feedrate will be expressed in respectively 
impacts per second (#/s) and Volt (V). To optimizse the parameter S , following 
10 quadratic criterion V{9) is proposed in function of prediction error s{tj S) (#/s): 



V{&) = £{0.5^^ (f, S)} e{t, S) = ;;(/) - y{t, S) (2) 

I r00491 I n case an optimal solution has to be found after N discrete input and 
15 output measurements (u^ yk) are carried out, V{S) can be approximated by 



vA^h\:jl^'^^"{k,3\ (3) 



20 [00501 F or each set of N measurements it is possible to find an optimal 
parameter 3 wherefor the quadratic criterion V{3) is minimised. Figure 2 gives 
the optimizsed exponential relation for five sections in one 330m harvest strip 
with a large variance in threshability due to different nitrogen applications and 
slope (a) fluctuations. In the first and last section, some erratic data points are 

25 lying near to the X and Y-axis respectively because of the dynamic step responses 
when the harvester enters and leaves the crop. 



[00511 R ecursive estimation of grain loss curves 
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As can be seen in Figure 2, a large variance in the model parameter can be 
detected in one single crop strip. Therefore, it is preferable to have an online 

11 



optimal estimation (in accordance with the invention) of the local separation 
behaviour. This site-specific information offers new possibilities for automatic 
tuning systems in a more direct way compared to the prior art online yield 
mapping systems, since it immediately relates feedrate u(t) with straw walker 
5 losses y(t) and provides an extra crop parameter i9 that may play an important role 
in the evaluation of new crop varieties. 

I [00521 Herein an online optimizsation procedure is derived based on a 
stochastic gradient method. The stochastic gradient method can be seen as a 
10 stochastic analog of the method of steepest descent for the minimisation of a 
deterministic function. In its general form, this method of steepest descent is 
represented by: 
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dx 



(4) 



where x^^^ denotes the t th iterate and y^^^ a chosen positive scalar. This approach 
may show useful , however, when the iterates are getting close to the minimum, 
this method is known to be fairly inefficient. 
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[00531 The so-called quasi-Newton methods yield distinctiy better results and 
use a modified search direction from the negative gradient direction: 



dx' 



V{x) 



dx 



V{x) 



(5) 



25 This iteration will provide convergence in one step to the minimum of V{x), if 
this fiinction is quadratic in x and /^'^ = 0.5 . 



[00541 W hen applied to the present optimizsation problem, the quasi-Newton 
optimizsation scheme can be transformed into following gradient scheme, which 
30 could be called a "stochastic Newton algorithm*': 
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V{3)^EH{3,e{t)) 



(6) 



3{t) = 3{t - 1) + ridv (^(/ - 1), ^' r Q(3{t - 1), ^0 (7) 

5 

where V" (.) denotes the approximate Hessian, -Q (.) is the gradient of H(x,e) with 
respect to x and indicates that the approximation depends on previous noise 
values = eft), eft-1), ... When this scheme is applied to the problem definition 
of equation (1) and (2), following algorithm is obtained for the model parameter 
10 &{t): 

3{t) = &{t - 1) + r{t)r-' {t}y(t, &{t - 3{t - 1)) (8) 

where i//{t,S)- ^^^^'^^ and scalar rft) corresponds with the instantaneous 
dS 

15 estimation of Hessian F"(i9). When an exponential process model according to 
function (1) has been chosen, ii/{f,9) will be equal to exp(i9.w(0).w(0 ■ 

[00551 T he accuracy of this approximation of the true Hessian plays a 
particularly important part when the recursive algorithm is operating close to the 
20 minimum. The natural approximation of this Hessian is to introduce the sample 
mean value: 



25 I [00561 H owever, the estimation of rft) in equation (9) puts as much attention 
on measurements temporally far firom t as on more recent measurements. 
Therefore, a weighted estimation of rft) usually yields better results: 



(0 = Z P{t> (k> ^{k - 1))+ sith ( 1 0) 
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where ro denotes the initial estimation of the Hessian function and weighting 
coefficients fi{t, k) and S{t) should be chosen such that 



t^fiM+S{t)^l V/, (11) 



[00571 A standard way to define the weighting coefficients is given by 
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It is easy to verify that this choice of parameters fulfils condition (11). Constant 
gain parameter corresponds to an exponential forgetting factor A© = ^ *" T^o • 

I r00581 A preferred but non-limiting practical realisation of this stochastic 
15 Newton algorithm is shown hereinafter: 



r (0)=ro; 
e(0)=eo; 
k-1; 

20 while (k<N); 

i|;{k)=u(k) *exp[u(k) *e(k-l) ] ; 
e (k)=y (k)-exp[u(k) *0(k-l) ]+l; 
r (k) =r (k-1) +Yo* [i^^ {k)-r(k-l) ] ; 
e(k)=e(k-l)+Yo*^(k)*e(k) /r (k) ; 

25 k++; 
end; 



100591 The algorithm is computationally cheap and can easily be implemented 
in machine software. Three parameters have to be determined before starting the 
30 algorithm. The preferred but non-limiting default values that are used in this study 
are added between brackets. 
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1. Gain sequence y{k) is set at a constant value (0.2). This parameter 
immediately determines the tracking capabilities of the algorithm but has also 
an important influence on the variance on the estimation of parameter &{k) . A 
5 high constant gain brings about fast tracking dynamics to follow fluctuations of 
S{k), but introduces large fluctuations around the true 3{k), even when the 
true parameter does not vary in time. It depends on the application purposes 
and sample frequency whether a high or low gain should be chosen. 

10 2. The initial estimation ro (500) of the Hessian function determines the 
confidence of having an accurate initial estimation of 3{k) A high value of ro 
implies a long term effect of the initial estimation &q on the following 
parameter estimations. 

15 3. Different strategies can be introduced to choose the initial parameter estimation 
i9(, (8). The last parameter value of a previous harvest run can be used or a 

parameter estimation of a near by crop strip when the tracking algorithm is 
connected to a positioning system. 

20 I r00601 R esults and discussion 

The algorithm that has been developed in the previous section can be used to track 
the relation between all combinations of process (e.g. separation) and feedrate 
sensors. Hereabove, the relation between feedrate and straw walker losses has 

25 been studied in more detail, since it may be used in automatic tuning systems that 
control the straw walker losses by adapting the feedrate. When the static non- 
linear characteristic between both signals can be tracked online, a ground speed 
control system can be realised that keeps the straw walker losses at a predefined 
level. When such knowledge about the instantaneous separation behaviour would 

30 not be available, more conservative controllers have to be designed, resulting in 
lower control performances. 
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[00611 In Figure 2, the signals of one 330 m harvest strip were divided into 
five distinct sections of about 66 m each. For each section, an optimal exponential 
I parameter 9 is calculated based on offline optimizsation techniques. Figure 3 
shows the results of the recursive algorithm for the same harvest strip. The results 
5 I of the previously described offline optimizsation procedure are also shown to 
illustrate the performance of the tracking algorithm. Roughly the same parameters 
are obtained, but now in an online, recursive way. The zone between 50 and 150 
m corresponds to data from an uphill section (+15%). Harvesting uphill reduces 
the separation capacity of conventional harvesters and as a result, larger 
10 coefficients 9 are obtained. Similar results are found for zones with smaller 
grain/straw ratios or higher moisture contents. 

[00621 S ite-Specific interpretation of estimated parameter sequence 

15 In case variations of parameter 9 depend on local field conditions, the same trend 
should be visible in adjacent strips. Figure 4 shows this type of parameter map for 
a set of 7 runs that are harvested uphill, parallel to the Y-axis. The correlation 
with the field slope or machine slope is evident, illustrating the value of this type 
of parametric maps for on-line tuning of harvesters, before variations in local 

20 processability are actually registered. This field was fertilised with different doses 
of Nitrogen in strips parallel to the Y-axis. Therefore, the harvester reacts in a 
slightly different way according to the Nitrogen gift and by this, to its position on 
the X-axis. A field with a constant Nitrogen application would give an even better 
correlation of the parameter and slope map. 

25 

I [00631 A lthough in the foregoing discussion the relationship between feedrate 
u(t) and straw walker losses y(t) is described, the method of the invention could 
I equally well be used for optimizsing parameters 9 in models for other crop 
handling processes at any of the following sections of a combine harvester: 
30 the separation section, e.g. the straw walker or the separator rotor; 

the sieves; 

the return flow system; or 
the cleaning section; 
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or indeed in sections of other kinds of harvester. 

I r00641 processability parameter can also be used for predicting the 

behaviour of the harvesting machine with respect to other effectiveness values 
5 such as the proportion of damaged useable crop parts, e.g. broken grain kernels, 
or a proportion of unwanted material in the useable crop parts, e.g. chaff and 
straw and particles in the clean grain. 

[0065] T he ability of the method of the invention to optimizse the parameter 
10 j9 in an online way offers numerous potential advantages, such as: 

the ability to devise a reliable vehicle control system that controls a 
performance variable of the harvesting machine and that relies on a feedrate set 
point calculated by an inverted form of Equation (1); 

the ability to devise processability-based expert systems, e.g. diagnostic 
15 systems that are self-executing in a combine harvester; warning systems 
manifested as feedrate or parameter "out of range" values; or information systems 
such as gradient-related data about likely crop losses. 

I £0066]__A further advantage of the methods of the invention is that the 
20 processability parameter 9 can if desired be assessed largely independently of the 
various sensing subsystems of a harvesting machine equipped for precision 

j farming. In other words the optimizsed threshability parameter 9 implicitly takes 
account of the crop, vehicle and field conditions that give rise to a particular value 
of S , without necessarily having to evaluate each individual cause. 

25 

I 100671 T he method described in full detail above can be used for optimizsing a 

single parameter used in a model for crop processing. However, it is readily 
I conceivable that analogous methods can be used for optimizing two or more 

stochastic parameters where such plurality of parameters is used for modelling an 
30 operation in a harvesting machine. Each parameter by itself, or each combination 

of parameters, can constitute a characteristic of the readiness with which the crop 

is processed. 
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CLAIMS 



What is claimed is: 



1 

2 1. A method of substantially continuously optimizsing a stochastic parameter 

3 3 that characterizses the instantaneously prevailing readiness with which crop is 

4 processed in a harvesting machine, including the step of recursively calculating 

5 I the optimizsed parameter value in accordance with the following algorithm: 
6 

7 m^fikt-'^u^tMt-m - (a) 

8 wherein: 

9 9{t) is the optimizsed stochastic parameter value at time t; and 

1 0 £{t^ 9{t)) is an error prediction function. 

1 
2 

3 2. A method according to claim 1, wherein G haracterised in that the algorithm 

4 (A) has the form: 

5 ^(0 = f[9{t - 1),..., 3it - e{t - n,u] 
1 

2 

3 3. A method according to Claim 1 or 2 , characterised w hcrein in that the 

4 algorithm (A) has the form: 
5 

6 &{t) = 3{t - ihrity-' {t)v^(t. &{t - i)>M(/ - 1)) 

7 

8 wherein: 

9 y(t) is a gain term; 

10 r(t) is a scalar approximation of a Hessian K"(»9) in which F is a 

1 1 quadratic error criterion; 
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(//{t,S) = ^^^^^ , in which y{tyS) is an estimation of a value indicative of 



13 the effectiveness of crop processing in said harvesting machine, said estimation 

14 being based on stochastic parameter 3 ; and 

15 €(t,3(t-l)) is the difference between the actual effectiveness value y(t) 

16 and the estimated value y{tyS) based on the previously optimizsed parameter 

17 .9(^-1). 

1 

2 

3 4. A method according to Claim 3, characterised w herein in that the 

4 algorithm (A) includes an estimation of r{t) that is weighted to reduce the 

influence, on the optimizsed parameter values i9, of past measurements. 

1 4t 
2 

3 5. A method according to Claim whcrein charaoterised in that : 

4 said stochastic parameter & is usable in a model for the relation between a 

5 value u(t) indicative of the feedrate of crop into the harvesting machine and a 

6 value y(t) indicative of the effectiveness of an operation processing said crop in 

7 said harvesting machine; and 
said value y{t, S) is an estimation value of the effectiveness obtained by 

the application of said model to the feedrate values u(t). 

1 
2 

3 6. A method according to Claim 5, charact e rised w herein in that said model 

4 comprises an exponential function. 



7. A method according to Claim 6, characterised wherein in that said model 
has the form: 

19 
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1 

2 

3 8. A method according to any of the Claims 54e-^, wherein c haracterised in 

4 that: 

5 said crop processing comprises separating useable crop parts from other 

6 plant matter; and 

7 said value y(t) is indicative of a flow of useable crop losses in a selected 

8 part (13/11) of the harvesting machine. 



9. A method according to any of the C laims 5-te-^, wherein€ 
that* 

5 said crop processing operation comprises separating useable crop parts 

6 from other plant matter; and 

7 said value y(t) is indicative of a flow of return crop in a selected part 
of the harvesting machine. 



1 

2 

3 10. A method of operating a harvesting machine comprising the steps of: 

4 (i)-substantially continuously optimizing a stochastic parameter & that 

5 characterizses the instantaneously prevailing readiness with which the harvesting 

6 machine processes crop: and 

7 (H)-substantially continuously adjusting a performance variable of the 

8 harvesting machine in dependence on the instantaneous, optimizsed value 3 of 

9 the-said parameter in order to optimizse the load of the harvesting machine so as 

10 to keep a value y(t) indicative of the effectiveness of said harvesting machine 



below a predetermined value. 



11. A method according to Claim 1 0, wherein o haracterised in that : 

20 



4 processing the crop comprises separating useable crop parts from other 

5 plant matter; 

6 I optimizsing the load of the harvesting machine comprises optimizsing the 

7 feedrate u(t) of crop into the harvesting machine; and 

8 the effectiveness value comprises losses y(t) of useable crop parts. 

1 

2 

3 12. A method according to Claim 10 or 11 , wherein c haracterised in that the 

4 step ff ^ of continuously optimizing a stochastic parameter S -includes carrying 

5 out the method steps of any of — the Claims 1 — te 9. 

2 

3 13. A method according to an)^ of the C laims 10 to 12 , wherein characterised 

4 in that t he step («)K)f adjusting a performance variable of the harvesting machine 

5 occurs in dependence on the output of an inverted form of a yield loss estimation 

6 function: 

7 i)(/,5) = exp(A^(0)-l. -(B) 

1 
2 

3 14. A method according to any of the C laims 10 to 13 , vyherein^ 

4 in that adjusting a performance variable comprises adjusting the travel speed of 

5 said harvesting machine or the actual cutting width of a header of said harvesting 

6 machine, 

2 

3 15. A method of mapping one or more field lots for variations in a stochastic 

4 I parameter & that characterizses the instantaneously prevailing readiness with 

5 which crop is processed in a harvesting machine, the method comprising the steps 

6 of: 
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(i)— operating a harvesting machine to harvest crop in a sa4d-field lot; 
(fi)"Simultaneously measuring the machine load and the machine 
effectiveness and determining the position of the machine in the field lot; 

(fH)-storing data indicative of the position of the harvesting machine at 

time t\ 

^+v)-using the measured machine load data u(t), and machine effectiveness 
data y(t) in an optimizsation of tbe-said parameter 9 ; and 

(¥)-mapping the optimizsed parameter values 9 obtained from the step of 
using the measured machine load data u(t) and machine effectiveness data v(t) in 
an optimization of said parameter 3 : <4V>-so as to produce a parameter map of the 

field lot. 



16. A method according to Claim 15, wherein c haract e ris e d in that the step of 
using the measured machine load data uft). and machine effectiveness data vfO in 
an optimization of said parameter S (4V)-includes carrying out an optimizsation 

according to any of th e Claims I4e4^. 

17. A method of operating a harvesting machine comprising the steps of: 
(i)— substantially continuously optimizsing a stochastic parameter 9 that 

characterizses the instantaneously prevailing readiness with which the harvesting 
machine separates useable crop parts from other plant matter: and 

(H)-sending a display signal, that is indicative of the instantaneous 

parameter value »9 , to a display device. 

18. A method according to Claim 17, wherein characterised in that the step (i) 
of optimizsing a stochastic parameter 9 includes carrying out the method of any 
of the Claims l4e-^. 
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19. A method according to Claim 17 . wherein or 18. charact e ris e d in that the 
display signal indicates an abnormal parameter value 3 . 

+^ 

20. A mM ethods according to Claim l a nv of the preceding cla i ms , wherein 
characteris e d in that s aid harvesting machine is a combine harvester and the crop 
is a grain-bearing plant. 

21. A mM ethods according to Claim 8 or 9 or to any other Claim referring 
thereto , wherein characteris e d in that the said selected part of the harvesting 
machine is selected from : 

the straw walkers-44^; 
the rotary separator; 
the sieves-444); 
the grain elevator; 
the return flow system-fl^; 
the cleaning section; or 
the axial threshing and separating rotor; 
of a combine harvester. 
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ABSTRACT 

METHODS OF OPTIMIZSING STOCHASTIC PROCESSING 
PARAMETERS IN CROP HARVESTING MACHINES 



The flow of crop material through a harvesting machine, such as a combine, can 
be described using a wide range of models, providing a relationship between the 
harvester load u(t), e.g. the mass flow at the inlet, and an effectiveness value y(t), 
e.g. the grain loss at the outlet. However, prior art models are not applicable to a 
wide range of harvesting conditions or are very complicated, requiring a multitude 
of inputs. The invention proposes to use a simple model, comprising a stochastic 
parameter & , which is continuously optimizsed to adjust the model to the latest 
prevailing working conditions. Such parameter may be considered to constitute a 
variable which characterizses the instantaneously prevailing readiness of the 
harvesting process. Such variable S can be used for establishing field maps 
showing the evolution of the harvesting operation itself. It can also be used in 
automatic routines which adjust the harvester speed in order to limit the grain 
losses. 
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